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This paper presents a recent study of a damaged generic transport model as part of a NASA research 
project to investigate adaptive control methods for stability recovery of damaged aircraft operating in off- 
nominal flight conditions under damage and or failures. Aerodynamic modeling of damage effects is performed 
using an aerodynamic code to assess changes in the stability and control derivatives of a generic transport 
aircraft Certain types of damage such as damage to one of the wings or horizontal stabilizers can cause 
the aircraft to become asymmetric, thus resulting in a coupling between the longitudinal and lateral motions. 
Flight dynamics for a general asymmetric aircraft is derived to account for changes in the center of gravity that 
can compromise the stability of the damaged aircraft. An iterative trim analysis for the translational motion 
is developed to refine the trim procedure by accounting for the effects of the control surface deflection. A 
hybrid direct-indirect neural network, adaptive flight control is proposed as an adaptive law for stabilizing the 
rotational motion of the damaged aircraft. The indirect adaptation is designed to estimate the plant dynamics 
of the damaged aircraft in conjunction with the direct adaptation that computes the control augmentation. 
Two approaches are presented: 1) an adaptive law derived from the Lyapunov stability theory to ensure that 
the signals are bounded, and 2) a recursive least-square method for parameter identification. A hardware-in- 
the-loop simulation is conducted and demonstrates the effectiveness of the direct neural network adaptive flight 
control in the stability recovery of the damaged aircraft. A preliminary simulation of the hybrid adaptive flight 
control has been performed and initial data have shown the effectiveness of the proposed hybrid approach. 
Future work will include further investigations and high-fidelity simulations of the proposed hybrid adaptive 
flight control approach. 


I. Introduction 

Aviation safety research concerns with many aspects of safe, reliable flight performance and operation of today’s 
modern aircraft to maintain safe air transportation for the traveling public. While air travel remains the safest mode of 
transportation, accidents do occur in rare occasions that serve to remind that much work is still remained to be done 
in aviation safety research. American Airlines Flight 587 illustrates the reality of hazards due to structural failures 
of airframe components that can cause a catastrophic loss of control. 1 Not all structural damages result in a loss 
of control. The World War II aviation history filled with many stories of aircraft coming back home safely despite 
suffering major structural damage to their airframes. Recently, the DHL incident involving an Airbus A300-B4 cargo 
aircraft in 2003 further illustrates the ability to maintain a controlled flight in the presence of structural damage and 
hydraulic loss. 2 

In damage events, significant portions of the aircraft’s aerodynamic lifting surfaces may become separated and as 
a result this may cause the aircraft’s carefully designed center of gravity (C.G.) to shift unexpectedly. The combined 
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loss of lift, mass change, and C.G. shift can manifest in an unstable, off-nominal flight condition resulting from the 
aircraft being out of trim that can adversely affect the ability for an existing flight control system to maintain the 
aircraft stability. In some other instances, aircraft’s damaged structures may suffer losses in structural rigidity and may 
develop elastic motions that can potentially interfere with an existing flight control system in an unpredictable manner. 
Moreover, the load carrying capacity of damaged structures may also become impaired and therefore can potentially 
result in excessive structural loading on critical lifting surfaces due to flight control inputs by the unaware pilot Thus, 
in a highly dynamic and difficult off-nominal flight environment with many uncertainties caused by damages to the 
aircraft, the inner-loop flight control must be able to cope with complex and uncertain aircraft responses that can 
greatly challenge an existing flight control system. 

Flight control of damaged aircraft in off-nominal flight conditions poses significant technical challenges in many 
areas of disciplines including aerodynamics, structural dynamics, flight dynamics and control, as well as human fac- 
tors. Thus, a comprehensive investigation from the aircraft integrated system perspective is needed to research and 
develop adaptive flight control technologies that can be used to retrofit conventional flight control systems in order 
to enable aircraft to achieve safe flight objectives. This comprehensive investigation would provide an integrated 
approach to damage effect physics-based modeling and simulation, safety-of-fiight assessment, flight control and re- 
covery, and adaptive system verification and validation. Damage effect physics-based modeling generates a knowledge 
base for understanding the behavior of damaged aircraft performance in the areas of flight mechanics, aerodynamics, 
and structural dynamics to address system interactions among various Sub-systems such as aircraft dynamics, airframe 
structures, engines, and flight control actuators. Using this knowledge base, flight mechanics of the damaged vehicle 
can be evaluated by flight simulation to assess classes of damage that can be recovered with different types of flight 
control effectors. On-board modeling provides state assessments of damaged aircraft in-flight that can be used to aide 
pilot’s decisions and control. Adaptive flight control is a critical technology that enables damaged aircraft to recover 
post-damage flight stability. Research in neural network adaptive flight control provides a possibility for developing 
an effective damage adaptive control strategy that can adapt the damaged aircraft to changes in the vehicle stability 
and control characteristics. 3 Emergency flight planning and post-damage landing technologies are also, investigated 
to aide pilots with an intelligent decision support system to identify a suitable landing site, flight path planning under 
a reduced flight envelope, and ultimately a safe landing execution strategy. 4 While adaptive flight control has been 
much researched, it has not been universally adopted in the aviation industry due to a number of software and stability 
issues that are inherent with any adaptive flight control system. Certification of these adaptive flight control systems 
is a maj or hurdle that needs to be overcome. Thus, research in adaptive system verification and validation is needed in 
order to develop stability and certification requirements for new adaptive flight control methods. 5 

This paper focuses on die flight mechanics and the adaptive flight control of damaged aircraft. The damage nature is 
primarily due to changes in the aerodynamic configuration of the vehicle brought about by various modes of damage 
that include airframe and flight control surfaces. . Some of these types' of damage can cause a rapid loss of vehicle 
stability and control resulting from a significant loss of lift capability and/or a loss of control power. A flight dynamics 
model for a damaged aircraft is developed to account for various damage effects including changes in aerodynamics, 
mass, inertias, and C.G. A trim analysis is presented to enable a rapid estimation of new trim states to maintain 
the aircraft flight conditions. Damage adaptive flight control methods are developed to enable stability recovery of the 
damaged aircraft. Research irithe adaptive reconfigurable flight control pro vides a method of risk mitigation for certain 
types of damage. Recent advances in neural network, direct adaptive flight control provide a foundation for much of 
this research. 7, 11,12 A hybrid direct-indirect adaptive control method is proposed to extend the current capability of the 
neural network adaptive flight control. The hybrid adaptive control includes an indirect adaptive law that performs an 
on-line estimation of plant dynamics of the damaged aircraftiThe stability of this indirect adaptive law is established 
by the Lyapunov stability theory. An alternative approach is also presented whereby a recursive least-square method 
is used for the parameter identification process. 

II. Damage Effect Modeling 

A twin-engine transport-class generic aircraft is chosen as a platform for the damage effect modeling . We will 
refer to this notational aircraft as a Generic Transport Model (GTM). Fig. 1 is an illustration of the GTM. Damage 
to an aircraft airframe and/or control surfaces can cause the aircraft to be out of trim, which consequently can lead to 
dynamic upsets of the aircraft flight states. Understanding the aircraft aerodynamic characteristics during a damage 
event is critical to developing flight control strategies for the stability recovery of a damaged aircraft In order to assess 
the damage effects on the GTM, aerodynamic modeling is performed to estimate the aerodynamic coefficients, and the 
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stability and control derivatives of the damaged GTM for various damage configurations under consideration. Damage 
aerodynamic 'characteristics will then be incorporated into a flight dynamic model of the damaged aircraft that will be 
used to develop adaptive flight control strategies. A controllability study can be performed using the flight dynamic 
model to determine which damage configurations are controllable and those that cannot be controlled. 

The damage effect aerodynamic modeling is performed using a vortex-lattice code developed at NASA Ames 
Research Center. 6 This computational fluid dynamics modeling is capable of rapidly computing the aerodynamic 
characteristics and control sensitivity of the damaged GTM due to various flight control surface inputs. The damage 
effects are modeled as partial losses of the left wing, left horizontal stabilizer, and vertical stabilizer, as shown in Fig. 
1. Wing loss represents one of the critical modes of damage that is a current focus of the research. 



Fig. 1 - Generic Transport Model 





Fig. 2 - Aerodynamic Coefficients Due to Wing Loss at a = 12° and / 3 = 0° 

Fig. 2 shows the damage effect due to wing loss for the damaged GTM. The effect of wing loss can be seen as a 
significant source of loss of lift capability of a damaged aircraft as the lift coefficient can be reduced by as much as 
25% for up to a 50% span loss of one of the wings. Changes in the pitching moment coefficient are also a result of 
the wing loss. Changing lift and pitch moment causes the aircraft to be out of trim that leads to the inability for the 
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flight control system to hold -altitude and flight path angle. Moreover, the aircraft lateral motion becomes a factor as 
a significant side force and yawing moment develop without any aileron or rudder input This lateral motion causes 
the longitudinal and lateral motions of the damaged aircraft to couple, resulting in changes in the angular rates. The 
stability of the damaged aircraft can be regained if sufficient control powers are still available to overcome the rolling 
and yawing moments as well as to retrim the aircraft in the pitch axis. 






Fig. 3 - Control Derivatives Due to Wing Loss at a = 12° and 0 = 0° 

For an ideal symmetric aircraft, the aileron deflection effects the roll control with insignificant contribution to the 
lift and pitching moment coefficients. The effect of wing loss causes the aileron deflection to induce a change in the lift 
coefficient as well as a change in the pitching moment coefficients as seen in Fig. 3(aHb). The abrupt changes in the 
lift control and pitch control derivatives are due to the complete loss of one of the ailerons for a wing loss that extends 
beyond 25% span. The consequence of this is that the damaged aircraft would exhibit a pitch-roll coupling when the 
ailerons are deflected asymmetrically. To maintain a trim state, the flight control must compensate for the unwanted 
pitch motion with the elevators. The situation is similar for the elevator control as the effect of wing loss introduces 
a change in the side force coefficient and a change in the yawing moment coefficient as seen in Fig. 3(c)-(d). Thus, 
a deflection of the elevators would result in a pitch-yaw coupling that must be compensated within the flight control 
system by adjusting the rudder control accordingly. Because of the asymmetry, the general motion of a damaged 
aircraft is coupled in all the three axes. As a result, any adaptive flight control strategy must be able to effectively 
handle this cross-coupled effect. 

III. Flight Dynamics of Asymmetric Aircraft 

The longitudinal motion of a symmetric aircraft is typically symmetric with respect to the aircraft fuselage ref- 
erence line. The lateral motion is uncoupled from the longitudinal motion owing to the aircraft symmetry. For a 
damaged aircraft, the symmetry may no longer be preserved depending on the nature of the damage such as wing 
damage. The as ymm etry of the damaged aircraft thus causes the longitudinal motion and lateral motion to couple 
together. Furthermore, the C.G. is shifted away from the x - z plane. The motion of an asymmetric damaged aircraft, 
therefore, must be understood in order to evaluate any flight control design. To this end, we consider an asymmetric 
aircraft with a C.G. offset from some reference location as shown in Fig. 4. The reference location is a fixed point 
located at the coordinate (xq , ?/o, xo) on the aircraft which may be taken as the original C.G. of the undamaged aircraft 
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in order to maintain the same coordinate reference frame. The C.G. of the damaged aircraft can move relative to this 
fixed reference point. 



Fig. 4 - C.G. Shift Relative to Reference Point O 

The damage effect resulting from a wing loss creates a larger C.G. shift in the pitch axis y than in the other two 
axes as shown in Fig. 5. This results in an additional rolling moment that the flight control must be able to compensate 
for using the available control surfaces in order to maintain the damaged aircraft in a trim state. 



Fig. 5 - C.G. Shift due to Wing Loss 


A. Linear Acceleration 

To understand the effect of the C.G. shift, the standard equations of motiomfor flight dynamics of a symmetric aircraft 
must be modified to allow for the asymmetry. Assuming a flat-earth model for a rigid body aircraft, the force vector 
in the body-fixed reference frame of the aircraft is 


F B = rn 


dv 

dt 




r dm ) ■— W 


a) 
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where W = mg —sin# cos#sin<£ cos#cos$ j is the gravitational force vector, u? = £ p q r j is the 
aircraft angular rate vector, r = r 4- Ar is the position vector of the reference location such that f is the position vector 

of the C.G. and Ar = £ Ax Ay A z j is the displacement vector from the C.G. to the reference location. 

The aircraft mass is assumed to undergo a change so that 

m ■== m* 4- Am (2) 


where m* is the original mass of the aircraft and Am < 0 is the mass change due to damages. 

Assuming that the change in the mass of the aircraft is instantaneous, the force vector then becomes 

„ dv dw dAr 

Fb — m— 4- m— xAr-f mw x — W (j) 

dt at at 

where is the speed of the C.G. relative to the reference location which is assumed to be small relative to v and 
therefore may be neglected. 

Transforming from the body-fixed reference frame to the inertial, reference frame yields 


F == Fb -f mw x v (4) 

Expanding Eq. (4) gives 

X = m (it + qAz — f Ay — rv 4- qw 4- gsin#) (5) 

Y = m (i) — pAz + tAx dru — pw — ycos^sin^) (6) 

Z = m(w -bp Ay — qAx — qu -f pu — ycos#cosd>) (7) 


The angular acceleration terms appearing in Eqs. (5)-(7) are a result of the C.G. shift. Thus, the linear acceleration 
of an asymmetric aircraft is coupled with its angular acceleration. 


B. Angular Acceleration 

We consider the angular momentum vector in the body-fixed reference frame 

Ha=/( r x (oo x r)] dm 4- J { r x v) dm (8) 

Expanding this expression yields 

— Iu> 4- mAr x v (9) 

where I is the mass moment of inertia matrix with respect to the reference frame at the reference location. 

The time rate of change in the angular momentum gives rise to the moment equation in the inertial reference frame 


ft, an b ^auj , dv ^ . . 

M == — (*wxli 5 = — h mAr x — 4- up x ixj 4- mw x (Ar x v) 

dt dt dt 

Expanding Eq. (10) results in the following moment equations 


GO) 


L = Ixxjp - Ix y q ~ Ixz r + Ixypr — ftpg + (ft ~ Iyy ) ? r + ft ft ~ ? 2 ) 

4- m (gu 4-Tty) A# 4- m (in — git) Ay -m (H ru) Az (1 1) 


M — I x yp 4- Iyyq I yz f “b fyzPQ ~~ ^x yq^ i (Jxx / zz) P T ~r 7 jc 2 (p 

— m (zo 4- pt>) Az 4- m (pu 4- no) Ay 4- m (w — n;) Az (12) 

XX = -/ xz p - I yz q + I zz r + I xz qr - I yz pr + (I yy - I xz )pq + I xy ft ~ P 2 ) 

-f m ( v — pw) Ax —m (u 4- gio) Ay 4- m (pu 4- qv) Az (13) 

Equations (11)-(13) indicate that the C.G. offset effectively creates additional moments on the aircraft. Cross 
coupling in both the linear and angular accelerations are present. Thus, the longitudinal and -lateral motions of the 
aircraft are generally coupled and the aileron or elevator commanded input therefore will affect the aircraft motion in 
both stability axes. 


6 of 23 


American Institute of Aeronautics and Astronautics 



C. Aerodynamic and Propulsive Forces and Moments 


Assuming that the engine thrust vector is aligned with the s-axis of the aircraft, then the forces and moments due to 


aerodynamics and the propulsion are 

X = 5 T T max + {Cl + ACx) QS sin a — (C* D + AC D )QS cosacos/3 (14) 

Y = {Cl + A Cy) QS- {Cl + A Co) QS sin 0 ( 15 ) 

Z = - (Cl + ACx) QS cos a- - (Cl + A C D ) QS sin a cos/? (16) 

L = (Cf + A Q) QSc (17) 

M - (C£ + AGn) QSc + $rT ma * (z e - z&) (IS) 

iV = (C£ + A Cn) QSc + <^A T^maxye + S T T max y 0 (19) 


where (x e , ±y ei z e ) are the centers of thrust and the subscript * denotes the force and moment coefficients for the 
undamaged aircraft evaluated at the reference location. 

We assume that the left and right engines produce the same amount of thrust with a combined maximum thrust 
equal to and are symmetrically positioned with respect to the aircraft fuselage reference line. Then St where 
0 < St < 1 is the throttle position corresponding to a desired total engine thrust, and Sat where — | < Sat < § 
is the throttle differential position difference that results in a desired engine differential thrust equal to the left engine 
thrust minus the right engine thrust. The incremental changes in these coefficients due to damages are defined as 

AC = A C 0 +AG a a+AC 0 l3 + AC s 8 (20) 

where AC = [ C L - C£ C D - C* D Cy - Cl Ci-C? C m - C* m C n - C* f, S = [ <S 0 <5 r 

is the flight control surface deflection vector, the subscripts a, /?, and S denote the derivatives, and the subscript 0 
denotes the coefficients at a = 0 and 0 — 0 . 


IV* Trim Analysis 

The aerodynamic forces on asymmetric aircraft include a non-zero side force component that is generally not 
experienced on symmetric aircraft. For a steady flight, the side force equation becomes 

mpcosflsin# 4- (Cf + ACV ) QS - {C* D + AC D )QSsin0 = 0 ( 21 ) 

The side force trim for the asymmetric aircraft can therefore be accomplished by trimming the aircraft at a non-zero 
bank angle <j> with zero sideslip angle 0, However, this would result in a limitation in the bank angle in coordinated 
turn maneuvers. Another side force trim approach is to trim the aircraft level with zero bank angle 0 but at a non-zero 
sideslip angle 0. In either case, the aircraft would have to be trimmed in both the longitudinal and lateral directions 
simultaneously by searching for the steady state solution of Eqs. .(14) to (16) with 6 = 0 or 0 == 0> The trim analysis 
thus computes the trim values for the angle of attack a , bank angle <p or sideslip angle 0, and engine throttle position 
St as functions of the aileron deflection <5 a , elevator deflection 5 e , and rudder deflection S T for a given aircraft Mach 
number and altitude. We assume that the engine thrusts will be symmetric at all times so that Sat = 0. 

If an undamaged symmetric aircraft has a mass m* and is flying wing-level, i*e,, d>* — 0, with zero control surface 
deflection at the original trim angle of attack a * , sideslip angle 0* — 0 , and trim throttle position 5^ corresponding to a 
lift coefficient ( 7 £, drag coefficient C% and side force coefficient Cy — 0. Then for small changes in the aircraft mass 
and aerodynamic coefficients, we can determine the incremental trim angle of attack, bank angle, and throttle position 
to maintain approximately the same trim airspeed V and flight path angle 7 * by taking small but finite differences of 
Eqs. (14) to (16) and setting them to zero, thus resulting in 

A$rT maa H-(ACx 4- Cl, a Aa 4- Cl,^A0 -I- Cl, 6$) QS sino7-s -{Cl 4- Cx. a Aa 4- Cl,@A0 -f Cl, 6$) QS cos a* Aa 
— (A Cp 4- Co,aAa 4- Cd,^A0 4- Cd,s^) QS cos a* 4* {Cd 4- Cd, aAa + Gd,p A0 4- Cd,8$) QS sina^Aa 

— mg cos ( 7 * 4 - a*) Aa — Amg sin ( 7 * 4 - a*) ~ 0 (22) 
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(A Gy 4 Cy^Aa 4 Cy,@A0 4 Cy t $8) QS — (Cp 4 Cp^Aa 4 C p $ A 0 4 C Yi6 S) QSA0 

4 mg [cos ( 7 * 4 a*) — sin ( 7 * 4 a*) Aa] A <j> = 0 (23) 


— (ACp + Ci^ctAa 4 Cl^A/3 4 Cl, 6 8) QS cos a* 4- (Cl -4 Cx, )Q! Aa 4 . Cl,p A0 4* Cl,s5) Q S' sin a* Aa 

— {AC d 4- Cu,a Aa 4 Cp t pA0 4 Cl, 6 8) QSsin a* — (Cp 4 Cp^Aa 4 Cp t p A0 4 Cd,$8) QS cosa*Aa 

- mg sin ( 7 * 4 a*) Act 4 Amgcos ( 7 * 4 a*) — 0 (24) 

To find the trim bank angle at zero sideslip angle, we set A0 = 0 in the equations above. Equation (24) then is a 
quadratic equation in terms of Aa whose solution can easily be computed as 



(25) 


with 

a = (Cpa sin a* -Cp )a cos a*) QS " (26) 

b — [(Cl 4 Cl,s8 - Cp,&) sin. a* - (Cp 4 Cp^8 4 Cl }C c) cos a *] QS - mg sin ( 7 * 4 a 5 ") (27) 

c — - [(A Cl 4 Cl,s8) cos a* 4 (A Cp 4 Cp^S) sin a*] QS 4 Amgcos ( 7 * 4 a*) (28) 

From Eq. (23), we now find the trim bank angle 

mg [cos ( 7 * 4 a*) - sm ( 7 * 4 ck*) Aa] 

Finally, the incremental trim throttle position can be solved directly from Eq. (22). 

Trimming the damaged aircraft with bank angle will result in a reduced bank angle limitation. This would poten- 
tially affect the aircraft’s turn capability. Moreover, the aircraft will not fly wing-level which would not be acceptable 
for a landing approach. Therefore, the damaged aircraft can be trimmed alternatively with the sideslip angle. This will 
enable the aircraft to maintain a level flight but the control authority of -the rudder control surface will be reduced since 
it has to compensate for the non- zero sideslip angle. To obtain the trim sideslip angle, we set A <p = 0 in the Eq. (23) 
and solve Eqs. (22) to (24) simultaneously. 

In examining Eq. (23) with Ad) — 0, it is noted that if the undamaged aircraft is in a cruise phase at a minimum 
drag, then the trim sideslip angle for the damaged aircraft could be large if the damage develops a significant side 
force. Typically, it is not advisable to fly the aircraft at a high sideslip angle because of the stability issue. Depending 
on the extent of damages, an effective trim approach may be one that uses a combination of the trim bank angle and 
sideslip angle. 

In cases where the rudder control power is insufficient due to damages, then the engine differential thrust throttle 
position 8a,t could be used to provide an additional control effector to trim the aircraft in yaw. Using engine differen- 
tial thrust for yaw control requires examining the issue associated with a slow engine response relative to the responses 
of typical flight control surfaces. While in theory the engine thrust can be used to trim the aircraft in yaw, often by 
the time the engine thrust is adjusted differentially to the correct trim value, the aircraft may have reached a different 
dynamic state due to the loss in airspeed and or altitude if the damage condition is severe enough to cause the aircraft 
performance to rapidly deteriorate in its flight envelope. Because of the time scale difference between traditional flight 
control surfaces and engines, engine actuator dynamics must be accounted for in the overall flight control strategy. 

In addition to using the engine differential thrust as a control effector, other flight control surfaces can be used in 
an overall control redundancy design strategy. This investigation would examine the control effectiveness of a various 
combinations of flight control surfaces. For example, wing spoilers can be used for roll control and the wing flap 
extension or deflection can be used for pitch control. Some of these control surfaces may have different time latency 
characteristics such as wing flaps versus spoilers. In the control and stability analysis, actuator dynamic model of slow 
systems should be included in the overall flight dynamic modeling. 

The trim analysis shows that upon damage, the damaged aircraft would have to be retrimmed with a new trim 
angle of attack a = a* 4 Aa, new bank angle <j> .= A3 or new sideslip angle 0 = A0 , and new throttle position 

— 8? 4 A5y- The trim a, 3 or 0, and St are all functions of the flight control surface deflection 8 as well as the 
aircraft damage configuration. In general, the stability and control derivatives needed to retrim the damaged aircraft 
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are not known. Thus, it is necessary that these parameters be identified in flight by a parameter identification process. 
Assuming that the effect of damage on the aerodynamics can be estimated, then a trim strategy is to initially retrim the 
damaged aircraft with zero control surface deflection by setting 6 — 0 in Eqs. (22) to (24). Then, using the inner-loop 
rate-command-attitude-hold (RCAH) control, the control surface deflection for the damaged aircraft can be obtained. 
This allows the trim values to be refined. Depending on the nature of damage, the trim refinement may be repeated 
until the damaged aircraft becomes completely trimmed. 

V. Damage Adaptive Flight Control 

Most conventional flight control systems utilize extensive gain-scheduling in order to achieve desired handling 
qualities. While this approach has proved to be very successful, the development process can be expensive and often 
results in aircraft specific implementations. Over the past several years, various adaptive control techniques have been 
investigated. 7 Damaged aircraft presents a challenge to the conventional flight control systems because the aircraft 
dynamics may deviate from its known dynamics substantially due to a significant degradation in the flight performance 
of the damaged aircraft. This makes it difficult for the conventional flight control systems to cope with changes in 
the stability and control of the damaged aircraft. Adaptive flight control provides a possibility for maintaining the 
stability of a damage aircraft by means of being able to quickly adapt to uncertain system dynamics. Research in 
adaptive control has spanned several decades, but challenges in obtaining robustness in the presence of unmodeled 
dynamics, parameter uncertainties, or disturbances as well as the issues with certification, verification and validation 
of adaptive flight control software prevent it from being implemented in flight control systems. 8 Adaptive control laws 
may be divided into direct and indirect approaches. Indirect adaptive control methods provide the ability to compute 
control parameters from on-line neural networks that estimate plant parameters. 9 Parameter identification techniques 
such as recursive least squares and neural networks have been used in indirect adaptive control methods. 10 In recent 
years, model-reference direct adaptive control using neural networks has been a topic of great research interests. 11 " 13 
Lyapunov stability, theory has been used to establish robustness of neural network adaptive control to ensure that 
adaption laws for neural network weight updates are asymptotically stable. 

In the current research, we adopt the work by Rysdyk and Calise 11 to develop a neural network adaptive con- 
trol with dynamic inversion for damaged aircraft. The adaptive flight control is able to provide consistent handling 
qualities without requiring extensive gain-scheduling or explicit system identification for a damaged aircraft. This 
particular architecture uses both pre-trained and on-line learning neural networks, and reference models to specify 
desired handling qualities. Pre-trained neural networks are used to provide estimates of aerodynamic stability and 
control characteristics required for model inversion. On-line learning neural networks are used to compensate for 
errors and adapt to changes in aircraft dynamics. As a result, consistent handling qualities may be achieved across 
flight conditions and for different damage configurations. An architecture of the neural network adaptive flight control 
is shown in Fig. 6. Furthermore, we will extend this architecture to include an indirect adaptive control element that 
provides an on-line estimation of the true plant dynamics. The estimation approach is provided by an adaptive law 
based on the Lyapunov stability analysis. In addition, we also consider a recursive least square method for the on-line 
estimation. 



Fig. 6 - Direct Neural Network Adaptive Flight Control Architecture 
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A. Linearized Plant Dynamics 

First, we need to arrive at a linear dynamics of the damaged aircraft for the feedback linearization control To maintain 
airspeed and altitude, the damaged aircraft has to be retri mm ed using the trim method above. The damaged aircraft 
stability must be recovered by the RCAH controller. This results in control surface deflections necessary to maintain a 
desired angular rate command. To design a linear RCAH controller, we want to eliminate the linear acceleration terms 
in Eqs. (1 1) to (13) corresponding to the uncompensated damaged aircraft dynamics of linear motion resulting from 
damages. Combining Eqs. (5) to (7) with Eqs. (11) to (13) yields 

fxxP ?xy4 — Jxzf “b fxypT ~ fxzPQ 4 (Izz lyy) Q r fyz 7 ^) 

4 rn (i qv 4 rw) Ax - mpvAy — mpwAz ~ (C* 4 ACi) QSc (30) 


fxyP 4 fyyQ I yz f 4 fyzPQ fxyQT “1“ (Jxx % zz) P T 4 Ixz — r^) 

- mquAx 4m(pu4 rw) Ay - mqwAz = ( C * 4 A Cm) QSc 4 S T T max (z e ~ z 0 ) (3 1) 

fxzp— fyzQ ~i“ fzz*r 4 IxzQT lyzpr 4 { fyy Ixx) Pq 4 Ixy P ) 

- mruAx - mrvAy 4 m (pu 4 qv) A z = (C* 4 A C n ) QSc 4 S T T max y o (32) 


where 

ACi = A Ci + Cy~- _ m (cos 0 cos ^ - cos 6 sin <b— 

c c QS \ c c 

An An n Az , /-> Ax , m 9 ( a j. Ax , - a Az 

AC m = ACm “ C x ~rr 4 C z —zr 77 ^ cos@cos4>—~ ~rsm0 — 

c c Q5 \ c c 

A/y a ^ ^ Ay Ax mg ( . , , Ax . „Ay 

ACW ••= A( 7 n 4 — — CV,-— — 777; ( cos 9 4 sm 0 -^r 

c c <35 \ c c 

where C%,C y , and 6 ^ are X, F, and Z force coefficients normalized to the djmamic pressure force <35. 
The linear dynamics of the damaged aircraft is computed by linearizing Eqs. (30) to (32) 

duj 


(r + AI) — = (£* + Afi) d> + (f| + Af 2 )cr + (g* + Ag) 5 


(33) 

(34) 

(35) 

(36) 


where A — [ A p A q A r | is the angular rate vector, cr - j" Acc A/? A$ A(Jy j is the trim parameter 
vector, and 

/'-’f* sy* 

°l,p U Lq °Z,r 


f* = QSc 
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Equation (36) is the angular acceleration equation of the asymmetric aircraft which can be written in a state-space 
form as 

c % = (Pi + AFi ) & 4 (F 2 + AF 2 ) <r 4 (G-f AG) 3 (37) 


where F x = F^f*, F 2 = 1*”%*, G - F^g*, AP X ■= I” 1 (ff 4 Aft) - F x , AF 2 == I" 1 (f| 4- Af 2 ) - F 2 , and 
AG - I” 1 (g* 4 Ag) - G. 

Under ideal situations, the plant dynamics of an undamaged aircraft is assumed to be known. However, for a 
damaged aircraft, the plant dynamics become uncertain as the stability and control derivative matrices AFi, AF 2 , and 
AG are usually unknown. Consequently, the flight control needs to be able to adapt to the uncertain plant dynamics 
of the damaged aircraft. The angular acceleration vector Co of the damaged aircraft may be written as the sum of an 
ideal angular acceleration vector of the undamaged aircraft and a differential angular acceleration vector Ac o due 
to damage as 

co == cbi 4 A Co (38) 


The ideal, undamaged aircraft plant dynamics can be written as 


u?i — Fi.cj 4 F 2 <r 4 G5 (39) 

where the stability and control matrices for the undamaged aircraft Fi, F 2 , and G are assumed to be known. 


B. Direct Neural Network Adaptive Control 

The goal of the adaptive flight control is to be able to fly the damaged aircraft whose handling characteristics is 
specified by a reference model. The control adaptation must be able to accommodate damages using the available 
flight control surfaces. A reference model is used to filter a rate command vector co c into a reference angular rate 
vector and a reference angular acceleration vector co m via a first-order model 

4 co n co m = CO n CO c (40) 

where co n = diag (io p , co qy w T ) is the frequency matrix. 

The reference frequency parameters must be chosen appropriately in order to obtain a good transient response that 
satisfies position and rate limits on the control surface deflection. For transport aircraft, typical values of the reference 
model frequencies u p , u. q , and co r are 3.5, 2.5, and 2.0, respectively. 3 In cases when the reference model is over- or 
under-specified, the parameters of the reference model must be adjusted. The tuning of the reference model parameters 
can be performed using an adaptive-critic approach to ensure that the flight control can track the reference model in 
order to achieve desired handling qualities. 16 

The reference model angular rate vector co m are compared with the actual angular rate output Co to form a tracking 
error signal co e — co m — co. A pseudo-feed back control vector u e is constructed using a proportional-integral (PI) 
feedback scheme to better handle errors detected from the roll rate, pitch rate, and yaw rate feedback. The error 
dynamics, defined by proportional and integral gains, must be fast enough to track the reference model, yet slow 
enough to not interfere with actuator dynamics. Hie issue with the integrator windup during a control saturation is 
addressed by a windup protection which limits the integrator at its current value when a control surface is commanded 
beyond its limit. The pseudo-control vector u e is computed as 


u e = Kpco e 4 Kj- f co e dr (41) 

Jo 
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In order to ensure low-gain error handling performance, the error dynamics is designed with natural frequencies 
that match the reference model frequencies in the roll, pitch, and yaw axes. A damping ratio is chosen with ( p — ( g = 
£ r = 1 / V2 . These frequencies and damping ratio are incorporated into the proportional and integral gains as 

K P “ diag (2( p u>p , 2^ q a> qy 2( r u) r ) (42) 

K / = diag (wp, u%) (43) 

A dynamic inversion is performed to obtain an estimated control surface deflection command 6 to achieve a desired 
angular acceleration vector u> d using the known plant dynamics of the undamaged aircraft from Eq. (39) as 

6 = G -1 (ch d -F a c5 - F 2 <r) (44) 


assuming that B is invertible. 

In order for the dynamic inversion control to track the reference model angular acceleration rate vector d? m , the 
desired angular acceleration vector Co d is set to be equal to 

• .Wd - 0m +■ u e - n ad (45) 

where u ad is an adaptive control augmentation designed to cancel out the dynamic inversion error, so that in an ideal 
setting, the desired angular acceleration rate ib d is equal to the reference model angular acceleration rate u> m as the 
tracking error goes to zero asymptotically. 

Because the true plant dynamics of the damaged aircraft is unknown and is different from the undamaged aircraft 
plant dynamics as can be seen from Eq. (37), a dynamic inversion will result from the control surface deflection 5. 
This error is equal to 

e = & — u> d = w — F — F 2 cr — GS (46) 

Comparing with Eq. (37), we see that the dynamic inversion error can also be expressed in terms of the unknown 
plant dynamics due to the damage effects 

e == Ad? = AF X £ + AF 2 <t 4- AGS (47) 

Substituting Eq. (45) into Eq. (46) results in 

£ = ~d>e - Vie T U ad (48) 

Combining Eq. (41) with Eq. (48) yields 

e = Ae + B (n ad ~ e) (49) 

where e~ j* uj e dr J and 



0 

1 

1 

, B = 

0 

A = 

j 

-K f 

-K P J 


I 


The adaptive control augmentation vector n ad is based on a neural network adaptation law by Rysdyk and Calise 1 1 
that guarantees boundedness of the tracking error and of the network weights using a single-hidden-layer sigma-pi 
neural network 

u od = W T /3(Ci > C 2> C 3 ) (50) 

where f3 is a vector of basis functions computed using a nested Kronecker product with Ci , C 2 , C 3 as inputs into the 
neural network consisting of control commands, sensor feedback, and bias terms. 

The network weights W are computed by an adaption law, which incorporates an adaptation gain T > 0 and an 
e-modification term (jl > 0 34 according to the update law 

W = -r (/3e T PB + n ||e r PB|| W) (51) 

where the matrix P solves the Lyapunov equation A r P + P T A = — Q for some positive-definite matrix Q and the 
norm is a Frobenius norm. 
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The e-modification term provides a robustness in the adaptation law. 14 The update law in Eq. (23) guarantees the 
stability of the network weights and the tracking error. The proof of this update law using the Lyapunov method is 
provided by Rysdyk and Calise. 11 Solving for the matrix P with Q — I, the update law can be rewritten as 

W = -r(/3V-f /i||V|| W) (52) 

where 

V = (I + K7 1 ) + ~ j'vldrKJ 1 (53) 

VI. Hybrid Direct-Indirect Adaptive Control Concept 

While the direct neural network adaptive law has been extensively research and has been used with good successes 
in a number of applications, the possibility of high gain control due to aggressive learning can be an issue. Aggressive 
learning is characterized by setting the learning rate F high enough so as to reduce the dynamic inversion error rapidly. 
This can potentially lead to a control augmentation command that may saturate the control authority. Moreover, high 
gain control may also excite unmodeled dynamics of the plant that can adversely affect the stability of the adaptive 
law. To address this issue, we are considering a modification to the present direct adaptive law to include an indirect 
adaptive law that provides an opportunity to perform an on-line estimation of the plant dynamics of the damaged 
aircraft explicitly. We call this approach as a hybrid direct-indirect adaptive control concept. The indirect adaptive 
law will provide an estimated plant d 3 mamic$ that will be used in the dynamic inversion. If successful, the control 
command will result in a smaller dynamic inversion error so that the learning of the direct adaptation neural network 
can be reduced. An architecture of the proposed hybrid adaptive control concept is shown in Fig. 7. ha the current 
study, we are developing some initial indirect adaptive laws for the on-line estimation of plant dynamics based on the 
Lyapunov stability theory and also the well-known recursive least-square method. Future research still remains ahead 
to rigorously investigate this proposed concept followed by high-fidelity simulations. 



Fig. 7 - Hybrid Direct-Indirect Neural Network Adaptive Flight Control Architecture 

A* Indirect Neural Network Adaptive Control 

We would like to estimate the unknown plant matrices using a linear-in-parameter neural network approach as 

AFi = W 2X (54) 

AF 2 = Wj/3 a (55) 

AG = Wj/3 5 . (56) 

where the hat symbol denotes the estimated plant matrices and ,3^, (3 a , j3 6 are some neural network architectures that 
may not be necessarily the same as /3 for the direct adaptation neural network. 

The error dynamics now can be expressed as 

e = Ae-f BW t {3 - BW-,^5 - BWj/3 CT <x - BWj 0 s S - BAe (57) 
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where Ae < e is the residual error from the estimation of the plant matrices. 

We now propose the following adaptive laws for the estimation of W U) W ff) and 


= I^wFPB 

( 58 ) 

= r ff /3 £r cre T PB 

( 59 ) 

W 8 == T s /3 Je r PB 

(60) 


where F& , F a , T$ > 0 are the adaptation gains. 

The proof is as follows: 

We let W = W # + W,W w == W£ + W„, W ff . - W* -f W^, and W* = WJ + W* where the asterisk 
symbol denotes the ideal weight matrices to cancel out the residual error Ae and the tilde symbol denotes the weight 
deviations. 

The ideal weight matrices are unknown but they may be assumed constant and bounded to stay within a A- 
neighborhood of the residual error Ae so that 


A.« sup W* J (3 - W* 1 - W* 1 fi a cr - f3 s 6 - Ae 

uj.<r,a 


We define the following Lyapunov function 
V — e T Pe 4- tr 


( W T W WjW CT WjW 5 

I . . -f ; 4 — h 


- a; . 


Ts 


where P > 0 and tr (A) denotes the trace of a matrix A. 

The time derivative of the Lyapunov function is computed as 


V — e r Pe 4- e f Pe 4- 2tr 


\ r 1 r * 

Substituting Eqs. (57) and (51) into the above equation yields 


( w T W , , W %Ws\ 


(61) 


(62) 


(63) 


V = -e r Qe + 2e r PB (w* T + W T j 0 - 2e r PB (w*J + w£) f3 u u> 

- 2e T PB (w* r + W j) /3 CT cr - 2e r PB (w| r + Wj) ^5 - 2e T PBAe 


+ 2tr 


w£w M w^w, , wJw/ 


-W T /3e T PB - M W r !|e r PB|| (w* + w) + + 

' * Fu Fa- 


Ts 


We note that tr (AB) — tr (B A) , so that 


e r PBW r /3 = tr (e r PEW^] = tr (w r /3e r PBj 

e T PBW^ u w = tr (e T PBW^ u «) = tr (w£/3 w £e r PB) 

e r PBW lP a a = tr (e r PBW 1(5 a F) = tr (w£/3 ff <re T PB) 

e T PBWj/3 s 5 = tr (e r PBWj '^) = tr (wJ / 3 5 5e T PB) 
Also, by completing the square, we have 


2tr [-/iW r j|e T PB|| (w* + Wj 


—2(i ||e T PB[| 



W* 

"IT 


(64) 

(65) 

( 66 ) 

(67) 

( 68 ) 

(69) 
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Since = W w , W CT - W CT , and Wj = W* Eq. (64) then becomes 


V < -e T Qe + 2e T PBA - 2p, e T PB ! 


W* 


w 


w* 


2tr 


W 




- /3 w ufo r PB ] + Wj ( ^ - /3 a ae T PB ) + Wj [ ^ - /3 5 <5e J PB 

4 cr 


W 5 
T s 


Since IIBII = 1, we establish that 


e T Qe < p( Q) ||e|J 2 
e r PBA < p (P) ||e|| ||A|| 


lle r PB|: 


w* 

2 

W* 

2 

< P (P) l|e|| 

2 


(70) 

(71) 

(72) 

(73) 


where p (Q) and p (P) are the spectral radii of Q and P. 

In order to guarantee that V < 0, we require that the trace operator be equal to zero, thus resulting in the adaptive 
laws in Eqs. (58) to (60). In addition, we also require that 




■2p(Q) 


(74) 


The time rate of change of the Lyapunov functionis then strictly negative and therefore it would guarantee that the 
signals are bounded. We note that e, u?, cr, 6 £ £<*> but e £ £ 2 since 


e T Q edt < p 


(Q) r 

Jo 


||ef dt<V(0)-V ( t -* oo) + 2p (P) / ||e|) ]jA|| dt 


f 

Jo 


-2 p P (P) r 

Jo 


W* 


-f W 


w* 


dt < oo (75) 


Utilizing Eq. (74), we have 


V(t 


oo) < V (0) - 2^(P) 



dt < 00 


(76) 


Thus, the value of V as t — > oo is bounded. Therefore, we establish that W -4 0, 


0 , W c 


0, 


and |jw,j|| — » 0 imply |je|| -4 0 as t — > 00 . This means that the adaptive laws will result in a convergence of the 

estimated AFi, AF 2 , and AG to their steady state values. In practice, the inputs <j> = \ ui T a T S must be 
sufficiently rich that contain enough frequencies to capture all the plant d}namics. In order for the on-line estimation 
to converge the their correct values, the inputs need to be a persistent excitation (PE) class of signals such that if there 
exist ao , oq, To > 0 then 9 


a ol < 


rt-bTo 

Jt 


4> (t) 4> T (t) dr < oql 


(77) 


We can also “robustify” the adaptive laws similar to Eq. (52) to better handle unmodeled dynamics and distur- 
bances by adding an e-modification term 14 to Eqs. (58) to (60) as 


w u = Fu {p a ve T PB - ^ ||e T PB|| W„) 

(78) 

W a = r a Ace^PB — p,^ ||e T PB|| W a ) 

(79) 

W s = Ts [p s 6e T PB - M ||e T PB|[ Wj) 

o 

OO 
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in which case the time rate of change of the Lyapunov function becomes 


V < -2 y. jje r PB]j 


W* 


+ w 


-2Mcl|e T PBjj 


w: 

2 




-2^||e T PB[| 

The effect of the e-modification is to increase the negative time rate of change of the Lyapunov function so that 
as long as the effects of unmodeled dynamics and or disturbances do not exceed the value of V, the adaptive signals 
should remain bounded. The e-modification thus makes the adaptive law robust to unmodeled dynamics so that the PE 
condition may not be needed. 15 


W* 


— 


lle r PB|| 


|W 




(81) 


B. Recursive Least-Square Parameter Identification 

While the indirect adaptive laws above provide a computational method for on-line estimation of the plant dynamics, 
it would be incomplete to not consider the well-known .least-square method which is equally robust in parameter 
identification process. If the dynamic inversion error is somehow can be estimated, then we should be able to apply & 
recursive least-square method to determine the weight matrices W^-, W a , and W $. Suppose the estimated dynamic 
inversion error can be written as 

e = & T $+Ae . ( 82 ) 

where = £ W£ Wj Wj j, 0 = j /3^u; (3 a cr (3 S 6 j , and As is the computational error in the 
estimated dynamic inversion error e, which may contain noise resulting from the on-line derivative computation of A 
since 

e = C> - F : w - F 2 <t - G S (83) 

where Co is the estimated angular acceleration which may be subject to computational errors. 

One method of computing Co is to use a backward finite-difference method 


0Ji= — Af- (84) 

to estimate Co at the i-th time step, but this method can result in a significant error if At is either too small or too large. 

Another approach is to collect n number of data points which will be used to generate an at least C l smooth curve 
in time using a cubic or B- spline method. This curve is then differentiated at their knots to find the estimated derivative 
values. In either case, the derivative computation will introduce an error source As. If the error is unbiased, i.e., it can 
be characterized as a white noise about the mean value, then the least- square method can be applied to estimate the 
plant dynamics. 

We consider a minimization of the following cost function 


J(w u , W ff , W,) = j* ||e - <3? r e| 2 dT 


(85) 


Our objective is to find recursive least- square adaptive laws for W w , , and To minimize the cost function, 

we compute the gradients with respects to the weight matrices, thus resulting in 


djJ 

<93? 


0i‘-* Te ) T 


dr — 0 


The recursive least-square formula using the gradient method is 


$ = RB (e T - 9 t §) 


( 86 ) 


(87) 


where 

R ■= -R00 r R 


( 88 ) 
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To show this, we see that from Eq. (86) 


f 6G T dT$> = f Oe T dr 

(89) 

Jo Jo 

Let t 

R -1 = f ee T dr > 0 
Jo 

(90) 

Then, differentiating Eqs. (89) and (90) results in 

R- 1 # + R- 1 # = 0e r 

(91) 

r - 1 = ee T 

(92) 

Also, we have 

RR" 1 ==!=* RR- 1 + RR -1 = 0 

(93) 


Substituting Eq. (92) into Eqs. (91) and (92) and solving for # and R yield the recursive least-square adaptive 
law. The matrix R is called the covariance matrix and the recursive least- square formula has a very similar form to the 
Kalman filter where Eq. (88) is a differential Riccati equation for a zero-order plant dynamics. We will show that the 
recursive weight update law is stable and results in bounded signals as follows: 

We let # — #'* ~f ^ with the hat and tilde symbols denoting ideal weights and weight deviations, respectively. 
Then, the error dynamics can be written as 

e < Ae + BW T J3 — + B A (94) 

We choose the following Lyapunov function 


L=y-ftx('$ r R- 1 #) (95) 

where V is the Lyapunov function for the direct neural network adaptive control and we have established that V < 0. 
The time rate of change of the Lyapunov function is computed as 

L = V + tr ^2# T R- 1 <§ + l> T R- 1 <§ ) (96) 

The weights <3? can be shown to converge to the ideal weights 3>* 9 so that 

# = ~H6G t & (97) 

Substituting Eq. (97) into Eq. (96) results in 

L — V - $ T 9e T § < 0 (98) 

Thus, the recursive least square weight update law 7 is stable. 

In practice, the recursive least-square method can be used to estimate the plant dynamics either continuously or 
discretely at every n data samples. Continuous time estimation requires solving the differential equations (87) and 
(88) at each time step. On the other hands, the discrete-time sampling estimation provides more flexibility in that the 
estimation can be executed after a specified number of data points have been collected. This would ensure that the 
signals contained in the sampled data are sufficiently rich to enable an accurate convergence. Another advantage of the 
recursive least-square method is that it provides an optimal noise filtering to minimize noise effects in the estimation 
of the plant dynamics. The discrete-time recursive least square formula is 

= ®k-i + ~ $fc-i] (99) 

R k “ Rfc-i - A _1 Rfc_i Ok (l + X~ 1 0^Hk-i0k^ 0fcRfc-i (100) 
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where k denotes the update cycle that repeats every n data samples, == k W J fc W^fc j is the weight 

matrix at k-t h update cycle, A is a forgetting factor that can be used to discount past data, and 
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(101) 


( 102 ) 


VII. Control Simulations 

The neural network adaptive flight control for damaged aircraft is evaluated in a medium-fidelity simulation test 
environment. The flight simulator is a fixed-motion simulator equipped with a pilot station, programmable displays, 
and a 120° field-of-view visual system as shown in Fig. 8. Pilot command inputs are received through a control stick, 
a rudder pedal, and a throttle quadrant. Flight control software includes a flight dynamics model of damaged aircraft 
as developed herein. Simulations are performed at a 30 Hz frequency. 




Fig. 8 - Flight Simulation Test Environment 

The damaged GTM is evaluated with various wing loss configurations. Fig. 9 shows the angular rates of the 
damaged GTM with and without the neural network adaptive flight control. The neural network control augmentation 
can be seen to quickly adapt to the changing dynamics of the damaged GTM. The roll, pitch, and yaw rates are quickly 
brought to zero to stabilize the damaged aircraft In contrast, without a neural network control augmentation, the 
aircraft rates are changing rapidly, particularly in the roll axis. A rapid increase in the pitch attitude can result in the 
damaged GTM reaching its stall angle of attack that would render the aircraft in a dangerous situation. 
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(a) 



Fig, 9 - Rate Control with and without Neural Network Adaptation 



t, sec 


Fig. 10 - Control Surface Deflection 

Fig. 10 shows the control surface deflections corresponding to the neural network control augmentation. The right 
aileron is commanded to move substantially to correct for a left turning rolling moment resulting from a left wing 
damage. A maximum aileron limit of 35° is nearly reached. Thus, it is possible that for certain damage scenarios, the 
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control augmentation will not be able to stabilize the damaged aircraft due to the control power limitation- In such 
situations, other types of control surfaces must be considered to provide additional control authorities for stabilization. 
Optimal control allocation approach must be incorporated into the neural network adaptive flight control to maximize 
the control effectiveness of all the available control authorities. 

To evaluate the hybrid adaptive flight control with the indirect adaptive law and the recursive least square method, 
a simulation was performed in MATLAB environment A damage configuration corresponding to a 30% loss of the 
left wing is selected. A step input pitch doublet is simulated. The tracking performance of the three control laws is 
compared in Fig. 11. 


(a) 





Fig. 1 1 - Pitch Doublet Tracking Performance 



Fig. 12 - Tracking Error Norm 
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As can be seen, the hybrid adaptive control with the indirect adaptive law is able to improve the tracking perfor- 
mance of the direct neural network adaptive control. The combined direct adaptive control with the recursive least 
square parameter identification actually outperforms both the direct and hybrid adaptive control approaches as the 
tracking error is significantly reduced as seen in Fig. 12. The control surface deflections to achieve this pitch maneu- 
ver are shown in Fig. 13 . The elevator deflection for this pitch maneuver is nearly saturated. The direct neural network 
adaptive control produces more overshoot than the hybrid adaptive control and the recursive least square approach. 
The left rolling moment is compensated by the right aileron input and the adverse yaw is compensated by a small 
rudder input For this simulation, the actuator dynamics is not included in the study. 


(a) 
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Fig. 13 - Control Surface Deflection 

VIII. Discussions 

The current research in the damage effect aerodynamic modeling has focused on single damage sites. Multiple 
damage sites can also exist in a damage event. The next step in damage effect aerodynamic modeling is to generate 
a predictive capability for multiple damage sites. An approach would be to conduct CFD modeling for representative 
multiple damage patterns. These modeling results are compared to single damage site models and a learning system 
can be developed to establish a surface response mapping between multiple and single damage sites. Using this 
damage surface response mapping, damage effects for any arbitrary damage pattern can be rapidly estimated. 

In the current adaptive flight control research, only traditional control effectors that include ailerons, elevators, and 
rudder are used. In severe damage situations, these control effectors may not be sufficient to stabilize and maintain 
good handling qualities of the damaged aircraft. Therefore, any adaptive flight control method must include a control 
allocation strategy that utilizes other potential control effectors that are otherwise not used in a conventional flight 
control system. These control effectors can include engine differential thrust for yaw control, wing spoilers for roll 
control, and wing flap extension or deflection for pitch control. Research in the areas of control redundancy design and 
reconfigurable control will be conducted to investigate optimal control allocation strategies for these control effectors. 
The issues of time-scale separation due to actuator dynamics are important for systems with different time latency 
such as engines and flaps and thus will be an area of adaptive flight control research. 17 

Damage effects can present a serious challenge to conventional flight control systems because the aircraft Sight 
dynamics may deviate from its nominal flight dynamics substantially as result of the degradation in the flight per- 
formance of the damaged aircraft. This makes it difficult for the conventional flight control systems to cope with 
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changes in the stability and control of the damaged aircraft. While neural network adaptive control offers a promise 
of being able to adapt to changes in flight dynamics of damaged aircraft, rigorous validation by simulations and flight 
testing will be pursued to explore areas of concern in the neural network adaptive flight control. One of the unresolved 
concerns is the learning characteristics of a neural network. If the dynamic inversion error is large due to a large dis- 
crepancy between the true and nominal plant dynamics used in the dynamic inversion control, die learning rate must 
be set sufficiently high in order for the neural network to reduce the error rapidly. As a consequence of the aggressive 
learning, the neural network tends to generate high gain control signals that may not be dynamically achievable. This 
potentially can cause a number of problems including control saturation, load constraints during flight being exceeded, 
excitation of unmodeled dynamics, and others. One potential solution is to introduce the proposed hybrid adaptive 
control that incorporates an explicit parameter identification based on an adaptive law derived from the Lyapunov 
stability method or a recursive least-square method to estimate the true plant dynamics, which would be used for the 
dynamic inversion control rather than the nominal plant dynamics. This approach potentially offers a way to reduce 
the dynamic inversion error that the neural network has to compensate for. 

Integrated flight dynamics modeling is another area research that addresses interactions among many types of 
physics problems during flight. An integrated flight dynamics model will be developed to include an aeroservoelas- 
ticity interaction model of a flexible-body vehicle dynamics with a propulsion model and its actuator dynamics. This 
integrated model will capture the combined effects of the 6-dof rigid body dynamics, structural dynamics of airframe, 
and propulsion model. Post-stall aerodynamics can be included in the aerodynamic coefficients and derivatives which 
can influence the flutter margin and the aerodynamic damping of the airframe. 

Structural interaction with a flight control system is critical to any flight control development. 18 ” 20 Elastic de- 
flection and mode shapes can adversely contribute to the vehicle stability and control, resulting in problems such as 
flutter, control reversal, structural frequency interaction within the flight control bandwidth, and others. Research 
in the area of aeroservoelasticity is very important for advancing the knowledge of damage adaptive flight control. 
Recent advances in fluid-structure interaction modeling using coupled computational fluid dynamics-finite element 
method provide a predictive capability for aeroservoelastic effects on the stability and control of damaged aircraft. 21 
New adaptive flight control methods will need to observe and obey structural load constraints imposed on a damaged 
airframe. The resulting adaptive flight control methods therefore would be more d 3 mamically achievable. Aeroser- 
voelastic frequency interaction with a safety-critical flight control system will be investigated in order to develop an 
integrated approach for dealing with potential issues with high frequency signals from elastic modes injecting into the 
frequency bandwidth of the rigid-body aircraft dynamics. Research in aeroservoelastic filtering and structural iden- 
tification for flight control will provide methods for assessing the elastic contribution of the airframe and developing 
adaptive flight control methods that can effectively filter out unwanted structural resonant modes within the flight 
control bandwidth. 


IX. Conclusions 

This paper has presented recent results on the modeling, control, and simulation of damaged aircraft as part of 
the aviation safety research at NASA. The damage effect aerodynamic modeling has been performed to provide an 
understanding of the control and stability of an asymmetric damaged aircraft. The effects of aerodynamic and control 
coupling in all the three stability axes are revealed from the modeling results. A 6-dof flight dynamics of symmetric 
aircraft is derived in order to account for the effect of the center of gravity shift resulting from the damage. An 
approach for trimming the damaged aircraft for the translational motion is presented. The trim procedure provides 
initial estimates of the trim values for the angle of attack, angle of sideslip, and engine thrust. The influence of the 
control surfaces on these trimmed values is then accounted for by adjusting the initial trim values with the control 
surface deflections obtained from the flight control. A hybrid direct-indirect neural network adaptive flight control 
concept has bees proposed to provide an opportunity to estimate plant dynamics in conjunction with the current direct 
adaptive control augmentation strategy. The on-line estimation of the plant dynamics is provided by an adaptive 
law derived from the Lyapunov stability theory and the recursive least-square method. The adaptive flight control is 
designed to track a reference model that specifies desired handling characteristics for a class of transport aircraft The 
feedback control augmentation uses a proportional and integral scheme to handle errors in the roll, pitch, and yaw rates . 
A control simulation of the direct adaptive control law is performed in a flight simulator to assess the stability recovery 
of a damaged generic transport model using the neural network adaptive flight control The results of the simulation 
show that the direct neural network control augmentation scheme is able to stabilize a damaged aircraft. In the near 
future, a control simulation of the hybrid adaptive control law will be conducted to investigate the potential benefits 

22 of 23 


American Institute of Aeronautics and Astronautics 



offered by this proposed scheme in reducing die possibility of high gain control in the present direct adaptive control 
strategy. Moreover, adaptive flight control research will advance the knowledge in the area of integrated flight control 
with propulsion and airframe effects in order to address interactions between vehicle dynamics, propulsion dynamics, 
and structural dynamics that may be present The assumption of rigid-body aircraft flight dynamics no longer holds 
true as the aircraft will have to be treated as an elastic body. This will give rise to challenges in developing adaptive 
flight control that can handle aeroservoelastic effects of damaged aircraft. 
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